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Abstract 

Transients  in  a  load  have  a  significant  impact  on  the  performance  and  durability  of  a  solid  oxide  fuel  cell  (SOFC)  system.  One  of  the  main 
reasons  is  that  the  fuel  utilization  changes  drastically  due  to  the  load  change.  Therefore,  in  order  to  guarantee  the  fuel  utilization  to  operate  within 
a  safe  range,  a  nonlinear  model  predictive  control  (MPC)  method  is  proposed  to  control  the  stack  terminal  voltage  as  a  proper  constant  in  this 
paper.  The  nonlinear  predictive  controller  is  based  on  an  improved  radial  basis  function  (RBF)  neural  network  identification  model.  During  the 
process  of  modeling,  the  genetic  algorithm  (GA)  is  used  to  optimize  the  parameters  of  RBF  neural  networks.  And  then  a  nonlinear  predictive 
control  algorithm  is  applied  to  track  the  voltage  of  the  SOFC.  Compared  with  the  constant  fuel  utilization  control  method,  the  simulation  results 
show  that  the  nonlinear  predictive  control  algorithm  based  on  the  GA-RBF  model  performs  much  better. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Compared  to  conventional  heat  engines,  steam  and  gas  tur¬ 
bine,  and  combined  cycles,  fuel  cells  are  considered  as  an 
efficient  electrical  power  generator.  Among  the  various  types 
of  fuel  cell,  solid  oxide  fuel  cell  (SOFC)  has  attracted  consid¬ 
erable  interest  as  it  offers  wide  application  ranges,  flexibility 
in  the  choice  of  fuel,  high  system  efficiency  and  possibility  of 
operation  with  an  internal  reformer  [1]. 

Transients  in  a  load  have  a  significant  impact  on  the  fife  of 
the  SOFC.  One  of  the  reasons  is  that  the  fuel  utilization  changes 
drastically  due  to  the  load  transient.  Fuel  utilization  uf  is  one 
of  the  most  important  operating  variables  of  the  SOFC  system. 
An  overused- fuel  condition  (wf  >  0.9)  could  lead  to  permanent 
damage  to  the  cells  due  to  fuel  starvation  and  an  underused- 
fuel  condition  (uf  <  0.7)  results  in  a  rapid  rise  of  the  cell  voltage 
[2],  For  protecting  the  SOFC  stack,  the  desired  range  of  fuel 
utilization  is  from  0.7  to  0.9. 

During  the  last  decades,  numerous  dynamic  models  have  been 
developed  to  study  the  effects  of  load  transients  on  the  perfor¬ 
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mance  and  durability  of  the  SOFC  system  [3-6].  These  models 
include  the  electrochemical,  thermal  and  mass  flow  aspects  of 
fuel  cell  operation.  Unfortunately,  very  few  articles  are  avail¬ 
able  on  the  development  of  effective  control  strategies  of  the 
fuel  utilization. 

The  terminal  voltage  and  the  fuel  utilization  cannot  be  kept 
constant  simultaneously  when  load  changes  [7].  So  there  are 
two  control  strategies,  which  can  guarantee  the  fuel  utilization  to 
operate  within  a  safe  range.  One  is  directly  controlling  the  input 
hydrogen  fuel  in  proportion  to  the  stack  current,  and  the  constant 
utilization  control  can  be  accomplished,  which  is  analyzed  in 
Ref.  [8].  The  other  one  is  to  maintain  a  constant  cell  voltage 
at  the  SOFC  terminal.  One  proper  value  for  cell  voltage  can 
guarantee  the  fuel  utilization  within  the  desired  safe  rang  when 
the  load  changes. 

In  this  paper,  a  nonlinear  model  predictive  control  (MPC) 
scheme  is  proposed  to  control  the  voltage  and  guarantee  the  fuel 
utilization  within  a  safe  range.  The  SOFC  system  is  an  uncertain 
nonlinear  system  and  its  structure  and  parameters  vary  with  the 
change  of  operating  point.  So  the  controller  should  be  robust  to 
uncertainty  and  meet  closed-loop  objectives  such  as  tracking, 
regulation  and  disturbance  attenuation. 

MPC  is  a  feedback  control  strategy  based  on  a  predictive 
model  and  receding  horizon  optimization.  The  important  advan¬ 
tage  of  MPC  comes  from  the  predictive  model,  which  allows  the 
design  of  controller  in  bigger  horizon  than  the  other  control 
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Nomenclature 

hi 

width  of  Gaussian  membership  function 

C'i 

center  of  Gaussian  membership  function 

CV 

conversion  factor 

E 

open-circuit  reversible  potential  (V) 

£° 

standard  reversible  cell  potential  (V) 

F 

Faraday’s  constant  (emol-1) 

I 

stack  current  (A) 

II 

limiting  current  (A) 

Aho 

valve  molar  constants  for  hydrogen  (mol  s-1  Pa) 

Ah2o 

valve  molar  constants  for  water  (mol  s-1  Pa) 

Ko2 

valve  molar  constants  for  oxygen  (mol  s-1  Pa) 

Kr 

constant 

n 

number  of  electrons  participating  in  the  reaction 

N 

number  of  cells  in  the  stack 

Pn2 

partial  pressure  of  hydrogen  (atm) 

Pn2o 

partial  pressure  of  water  (atm) 

Po2 

partial  pressure  of  oxygen  (atm) 

P 

cell  pressure  (atm) 

P  Ac 

rated  power  (kW) 

<7f 

nature  gas  flow  rate(mol  s-1 ) 

4h2 

input  hydrogen  flow  (mol  s-1) 

*?h2 

output  hydrogen  flow  (mol  s-1) 

<fn2 

hydrogen  flow  that  reacts  (mol  s-1) 

qf1 

input  flow  (mol  s-1) 

(fi 

output  flow  (mol  s-1) 

flow  that  reacts  (mol  s-1) 

«o2 

input  oxygen  flow  (mol  s-1) 

m-o 

hydrogen-oxygen  flow  ratio 

R 

gas  constant  (J  mol-1  K) 

^Ohmic 

Ohmic  resistance  (fi) 

T 

cell  temperature  (K) 

Ui 

membership  function 

Uf 

fuel  utilization 

U 

output  voltage  (V) 

Ur 

rated  voltage  (V) 

V 

compartment  volume  (m3) 

Greek  symbols 

a 

electron  transfer  coefficient  at  the  electrode 

P 

Tafel  slope 

O'dCt 

activation  losses  (V) 

V  cone 

concentration  loss  (V) 

^Ohmic 

Ohmic  losses  (V) 

m2 

response  time  for  hydrogen  flow  (s) 

TH2o 

response  time  for  water  flow  (s) 

ro2 

response  time  for  oxygen  flow  (s) 

Tl 

reformer  time  constants  (s) 

d 

Tafel  constant 

algorithms  without  the  form  of  predictive  model  [9].  There¬ 
fore,  the  model  is  the  key  to  determine  the  control  quality. 
The  traditional  approach  usually  is  based  on  approximate  lin¬ 
earization  theory,  which  is  often  difficult  to  identify  an  accurate 


mathematical  model  of  the  system  and  imposes  serious  restric¬ 
tions  on  the  structure  of  nonlinear  systems  [10].  Moreover, 
the  robustness  of  the  closed-loop  system  cannot  be  guaran¬ 
teed,  especially  when  the  parameters  of  plant  are  uncertain  or 
there  is  noise  or  disturbance  in  the  process  [11].  Therefore  spe¬ 
cific  nonlinear  modeling  approaches  might  be  required.  Neural 
networks  are  considered  as  an  attractive  structure  to  establish 
the  mathematical  relationship  of  the  dynamic  system  based 
on  the  input-output  data.  A  RBF  neural  network  is  a  feed¬ 
forward  neural  network  with  one  hidden  layer  and  can  uniformly 
approximate  any  continuous  function  to  a  prospected  accuracy 
[12],  However,  a  key  problem  by  using  the  RBF  neural  net¬ 
work  approach  is  about  how  to  choose  the  optimum  initial 
values  of  the  following  three  parameters:  the  output  weights, 
the  centers  and  widths  of  the  hidden  unit.  If  they  are  not 
appropriately  chosen,  the  RBF  neural  network  may  degrade 
validity  and  accuracy  of  modeling  [13].  To  assure  the  optimal 
performance  of  the  RBF  neural  network  approach  for  SOFC 
modeling,  we  consider  applying  a  genetic  algorithm  (GA)  to 
optimize  the  RBF  neural  network  parameters  in  this  study. 
Genetic  algorithms  are  a  kind  of  self-adaptive  global  search¬ 
ing  optimization  algorithm  based  on  the  mechanics  of  natural 
selection  and  natural  genetic  [14],  Different  from  conventional 
optimization  algorithms,  genetic  algorithms  are  based  on  pop¬ 
ulation,  in  which  each  individual  is  evolved  parallel,  and  the 
ultimate  result  is  included  in  the  last  population.  In  this  paper, 
a  physical  model  of  a  100  kW  SOFC  system  is  used  to  generate 
the  data  required  for  the  training  and  predicting  of  the  GA-RBF 
model. 

This  paper  is  organized  into  six  sections.  Section  2  briefly 
describes  the  SOFC  system.  A  dynamic  SOFC  mathematical 
model  is  built  in  Section  3.  In  Section  4  a  nonlinear  predictive 
control  algorithm  based  on  a  GA-BRF  model  is  explained  in 
detail.  In  addition  the  constant  fuel  utilization  method  is  also  ana¬ 
lyzed.  Section  5  shows  the  results  obtained  with  the  two  control 
algorithms.  The  last  section  concludes  the  paper  by  analyzing 
the  obtained  results. 


2.  SOFC  system  description 

The  SOFC  system  includes  a  fuel  processing  unit  or  the 
reformer  and  a  fuel  cell  stack.  Hydrogen  is  a  main  fuel  for  most 
type  of  fuel  cells.  Nevertheless,  other  fuels  such  as  methane, 
methanol,  ethanol,  gasoline  and  oil  derivatives  can  also  be  used 
when  a  reformer  is  included  in  a  fuel  cell  system  for  converting 
the  fuel  to  hydrogen. 

The  basic  components  of  the  SOFC  are  anode,  cathode 
and  two  ceramic  electrodes.  In  the  fuel  cell,  fuel  is  sup¬ 
plied  to  the  anode  and  air  is  supplied  to  the  cathode.  At  the 
cathode-electrolyte  interface,  oxygen  molecules  accept  elec¬ 
trons  coming  from  the  external  circuit  to  form  oxide  ions.  The 
electrolyte  layer  allows  only  oxide  ions  to  pass  through  and  at 
the  anode-electrolyte  interface,  hydrogen  molecules  present  in 
the  fuel  react  with  oxide  ions  to  form  steam  and  electrons  get 
released.  These  electrons  pass  through  the  external  circuit  and 
reach  the  cathode-electrolyte  layer,  and  thus  the  circuit  is  closed 


234 


X.-J.  Wu  et  al.  /Journal  of  Power  Sources  179  (2008)  232-239 


[15].  The  electrochemical  reactions  are  given  as  follows: 


Anode  side  : 

H2  +  02-^  H20  +  2e“ 

(1) 

Cathode  side  : 

l/202  +  2e_  — >  O2- 

(2) 

Balance  reaction 

:  H2  +  1  /202  — ►  H20 

(3) 

A  single  cell  produces  an  open-circuit  voltage  of  approxi¬ 
mately  1  V,  so  fuel  cells  have  to  be  connected  together  in  a  series 
arrangement  to  form  a  stack. 

Fuel  utilization  is  one  of  the  most  important  operating  vari¬ 
ables  that  may  affect  the  performance  of  fuel  cell.  It  is  defined 


_  4h2  ~  <?h2  =  <4  =  NI 
<7h2  <7h2  2  Fq%2 

3.  SOFC  system  dynamic  model 

3.1.  Fuel  processing  unit  model 


(4) 


In  Ref.  [16],  the  author  introduced  a  simple  model  of  a 
reformer  that  generates  hydrogen  through  reforming  methane. 
The  model  is  a  second-order  transfer  function.  Ref.  [7]  shows 
that  fuel  process  approximated  by  first-order  transfer  functions 
was  suitable  for  simulation  point  of  view.  So  the  mathematical 
form  of  the  model  can  be  written  as  follows: 


<  =  cv 

<7f  TlS+1 


(5) 


When  these  reactants  are  fed  into  the  SOFC  stack  where  the 
reaction  described  by  Eq.  (3)  occurs,  from  which  it  is  seen  that 
full  reaction  ratio  between  hydrogen  and  oxygen  is  2  to  1 .  In 
order  to  allow  for  oxygen  to  completely  react  with  hydrogen  and 
maintain  the  pressure  difference  between  the  anode  and  the  cath¬ 
ode  below  a  certain  threshold  value,  excess  oxygen  is  provided. 
And  oxygen  flow  rate  is  determined  using  the  hydrogen-oxygen 
flow  ratio  rn-o- 


3.2.  SOFC  stack  model 


3.2.1.  Electrochemical  model 

The  electrochemical  oxidation  of  CO  is  not  considered  in  our 
study,  because  data  in  the  literature  show  that  the  electrochemical 
oxidation  of  hydrogen  is  faster  than  that  of  CO  [17].  Hence,  the 
Nemst  voltage  E  is  mainly  dependent  on  the  electrochemical 
reaction  of  H2. 

In  terms  of  input  and  output  flow  rates  and  exit  molarity,  the 
zth  component  material  balance  for  the  fuel  cell  stack  can  be 
written  [18,19]: 


d Pi. 
d  t 


K1  j 

-y(9i  ~9 i  ~Qi) 


(6) 


According  to  the  basic  electrochemical  relationships,  the 
mole  flow  that  reacts  q]  can  be  calculated  as 


Replacing  the  mole  flow  that  reacts  by  Eq.  (7),  we  can  get 
the  following  expression  of  the  hydrogen  partial  pressure 
d  RT  in  n 

^Ph2  =  —  (q%2-q0H2-2KTD  (8) 

For  orifice  that  is  choked,  it  could  be  considered  that  the  molar 
flow  of  any  gas  through  the  valve  is  proportional  to  its  partial 
pressure  inside  the  channel,  according  to  the  expression 


<?h2 

Ph2 


Kh2 


(9) 


Replacing  the  output  flow  by  Eq.  (9),  applying  the  Laplace 
transformation  to  Eq.  (8)  and  isolating  the  hydrogen  partial  pres¬ 
sure,  yields  the  following  expressions: 


m(*)  =  -  2Ktn  (io) 

1  +  th2s 
where  th2  =  7^ 

Similarly  component  balances  for  O2  and  H2O  lead  to  the 
following  set  of  equations: 


Po2(s)  - 


1/^02 
1  +  Tq  2S 


(4  -  KrI) 


(11) 


Fh2oCs)  =  27fr/ 


1/^h2o 
1  +  m2o  s 


(12) 


3.2.2.  Operating  cell  voltage 

The  actual  cell  potential  is  decreased  from  its  ideal  potential 
because  of  several  types  of  irreversible  losses,  such  as  activation, 
concentration  and  Ohmic  losses.  The  equation  for  polarization 
is  shown  as  follows: 

7 act  =  a  +  £  log  /  (13) 

RT  f  I\ 

7conc  =  2F  1°  (l  -  7l  J  (14) 

70hmic  —  I  R  (15) 

Applying  Nemst’s  equation  and  Ohm’s  law  (taking  into 
account  Ohmic,  concentration,  and  activation  losses),  the  stack 
output  voltage  is  represented  as  follows: 

Vdc  =  E  —  7  Ohmic  —  7conc  —  7act  (16) 

where  £  =  «  (^  +  ff  In  ^ 

The  dynamic  physical  model  replaces  the  real  SOFC  stack  to 
generate  the  simulation  data  required  for  the  identification  of  the 
GA-RBF  model.  The  parameters  of  the  SOFC  system  are  given 
in  Table  1 .  The  data  sources  blocks  developed  in  MATLAB  is 
shown  in  Fig.  1 . 


4.  SOFC  control 


4.1.  Constant  voltage  control 


*-"-2  KJ 

'll  n.w  r 


(V) 


The  basic  frame  of  the  proposed  model  predictive  controller 
for  the  SOFC  system  is  shown  in  Fig.  2,  where  yr  is  the  refer- 
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Fig.  1.  SOFC  dynamic  model. 


Table  1 

The  blocks  developed  in  MATLAB 


Parameter 

Unit 

Value 

T 

K 

1273 

F 

Cmol-1 

96485 

R 

J(molK)-1 

8.314 

N 

- 

384 

£° 

V 

1 

Kr 

mol(s  A)-1 

0.996  x  10“3 

Kr2 

mol  (s  Pa)-1 

8.32  x  10“6 

Kh2o 

mol  (s  Pa)-1 

2.77  x  10“6 

Ko2 

mol  (s  Pa)-1 

2.49  x  10“5 

*02 

s 

2.91 

*h2 

s 

26.1 

*h2o 

s 

78.3 

Uphmic 

Q 

0.126 

«h2 

mols-1 

Variable 

P 

- 

0.11 

9 

- 

0.05 

It 

A 

800 

n 

- 

2 

U 

V 

Variable 

I 

A 

Variable 

* 1 

s 

5 

Ot-O 

- 

1.145 

CV 

- 

1 

UT 

V 

330 

Pdc 

kW 

100 

Ufs 

- 

0.8 

ence  track  of  the  control  system,  y  and  vr  is  the  SOFC  dynamic 
model  output  and  the  GA-RBF  predictive  model  output,  respec¬ 
tively. 

A  nonlinear  model  that  the  controller  will  use  for  the  opti¬ 
mization  is  the  first  step  in  designing  a  nonlinear  MPC  system. 
This  model  should  be  as  accurate  as  possible,  while  being 
simple  enough  to  allow  for  repeated  calculations  during  the 
optimization.  If  the  dynamic  model  developed  in  Section  3 
is  applied  in  the  predictive  control  scheme,  it  will  consume 
much  time  to  obtain  the  solutions.  For  this  reason,  a  nonlin¬ 
ear  offline  voltage  model  of  SOFC  is  built  by  a  GA-RBF  neural 
network. 


4.1.1.  Predictive  model  based  on  GA-RBF  neural  network 
Suppose  the  operating  temperature  and  pressure  of  the  SOFC 
are  kept  constant  in  this  paper.  The  oxygen  flow  is  expressed  by 
using  the  hydrogen-oxygen  flow  ratio,  so  the  terminal  voltage 
U  is  mainly  influenced  by  the  inlet  hydrogen  flow  and  cur¬ 
rent  I.  The  following  NARMAX  model  is  used  to  describe  the 
controlled  voltage  system 

y(k)  =  f[y(k  -  1),  y(k  -2 ) . y(k  -  ny),  u(k  -  1), 

u(k  —  2), ... ,  u(k  —  nu)\  (17) 


where  y  is  the  SOFC  terminal  voltage,  u  is  the  hydrogen  flow 
rate  and  current,  ny  and  nu  are  the  lags  of  the  output  and  input, 
respectively,  mdfi-)  is  a  nonlinear  function.  In  this  section  we 
adopt  a  GA-RBF  neural  network  to  identify  the  nonlinear  func¬ 
tion^-).  The  structure  of  the  RBF  neural  network  is  shown  in 
Fig.  3. 

The  output  of  hidden  layer  is 


Ui  =  exp 


(x  -  Ci)T(x  -  a)  1 

2  bj  \ 


(i  =  l, 2, ...,*) 


(18) 


where  x  =  (/,  )T ,  Ci  =  (c,i,  cq)t  is  the  center  of  the  ith  RBF 

hidden  unit,  and  £,■  is  the  width  of  the  ith  RBF  hidden  unit. 


Fig.  2.  Block  diagram  of  the  proposed  control  system. 
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The  output  is 

y R  =  (19) 

!= l 

where  yR  =  U\w-  [uq,  W2,. .  .w9]r  is  the  neural  network  weight 
that  connects  output  >’r  and  neuron  i  in  the  hidden  layer. 

In  our  case,  it  is  done  by  minimizing  the  following  quadratic 
cost  function  of  the  output  error.  This  error  is  calculated  by 
comparing  the  output  value  of  the  network  and  the  desired  output 
value. 

1  M 

J  =  2Y^[yR(k)  -  y{k)f  (20) 

k=  1 

where  yR  is  the  output  voltage  calculated  from  neural  network; 
y  the  output  voltage  of  the  dynamic  physical  model  developed 
in  Section  3;  M  the  number  of  training  data.  A  gradient  descent 
algorithm  is  adopted  to  minimize  /.  And  the  genetic  algorithm  is 
used  to  obtain  the  optimum  initial  values  of  the  following  three 
parameters:  the  output  weight  wu  the  centers  c,-  and  widths  hi. 

The  predictive  output  of  p  step  can  be  obtained  by  the  GA- 
RBF  model  as  [20]: 

yR(k  +  p)  =  fR[y(k  +  p-l),y(k  +  p-  2), ..., 

y(k  +  p  -  riy),  u(k  +  p  —  1),  u(k  +  p-  2), ..., 
u(k  +  p  -  «„)]  (21) 

And  the  error  at  instance  k  is  defined  as 
e(k)  =  y(k )  -  yR(k)  (22) 

Therefore  the  predictive  output  of  feedback  system  is 
y(k  +  j )  =  yR(k  +  j )  +  e(k)  (j  =  1 , 2 . p)  (23) 


of  output  voltage  are  introduced  to  avoid  excessive  movement 
of  the  control  input,  which  are  defined  as 

yr(k  +  j )  =  cjy(k)  +  (1  -  cj)r  (1  <  j  <  p,  0  <  c  <  1) 

(24) 

where  yr  and  y  are  the  reference  trajectories  and  predictive  out¬ 
put,  respectively,  r  is  the  system  set  value.  The  optimization 
problem  for  the  model  predictive  controller  is  the  minimization 
of  the  sum  of  squared  errors  between  the  referenced  trajectory  yr 
and  the  predictive  output  y  with  an  additional  penalty  imposed 
on  rapid  changes  in  the  manipulated  variables.  Define  objective 
function  as 

m  =  X>R(fc  +  j)  -  yr(k  +  j)]2  +  jhku2(k  +  1  —  1) 

7=1  7=1 

u{k  +  i  —  l)eU*  (25) 

where  n  and  m  are  the  predictive  horizon  and  control  horizon, 
respectively.  U*  =  [wmin,  Wmax]  is  the  control  range.  Golden  mean 
method  is  used  to  obtain  the  optimal  output  of  controller  u*(k), 
the  optimization  process  is  presented  as  follows  [20]: 

•  Step  1:  Define  the  initial  searching  area  [ari,  ]  =  [um;n, 
Wmax];  the  initial  number  of  steps  k=  1;  and  define  £, 
it  is  positive  and  adequately  small;  the  initial  searching 
points  are  given  as:  A.  i  =ci\  +  (1  -0.618)(/li  -ai),Aj  =  a\  + 
0.618(/3|  —  ai).  And  then,  compute  the  values  of  objective 
function  J(/.])  and  J(k\). 

•  Step2:  If  fit  ~  otk  <  L  then  u(k)  =  (Jik  +  ak)/2;  Otherwise,  if 
J(Xk)  >  J(X'k),  go  to  step  3;  if  J(Xk)  >  J(X'k),  go  to  step  4. 

•  Step  3:  ak+i  =  Xk,  fik+\  =  fik,  Xk+\  =  k'k,  then  X'k+1  =  ak+i  + 
0.61 8(^+1  -  ak+\)\  and  then,  compute  the  value  of  objective 
function  J(X'k+]),  go  to  Step  5. 

•  Step  4:  ak+i  =  ak,  pk+i  -  l'k,  X'k+1  —  lk,  then 
Xk+i  =ayt+i+(l  -0.618)(jSyt+i  -ak+i)-,  and  then,  com¬ 
pute  the  value  of  objective  function  J(Xk+] );  go  to  Step 
5. 

•  Step  5:  k=K+  1,  return  to  Step  2. 

After  iterations  for  several  times,  the  optimal  control  moves 
u*(k )  can  be  obtained. 

4.2.  Constant  fuel  utilization  control 


where  p  is  the  predictive  horizon.  The  method  of  the  offline 
model  is  simple  and  convenient  only  if  the  experimental  data  or 
experience  of  the  input  and  output  variable  is  obtainable.  It  is 
not  necessary  to  decide  the  coefficient  of  the  SOFC  material  and 
structure  of  the  mechanism  model.  Therefore  it  can  be  used  in 
the  control  system. 

4.1.2.  Optimization  algorithm 

At  each  sample  time  k,  a  series  of  future  output  values  y(k  +  j ) 
is  calculated  through  GA-RBF  predictive  model  and  it  is  com¬ 
pared  to  the  reference  output  value  yv(k).  Referenced  trajectories 


According  to  Eq.  (4),  the  operation  of  the  SOFC  stack  with  a 
fuel  input  proportional  to  the  stack  current  results  in  a  constant 
utilization  factor  in  the  steady-state.  Thus,  the  SOFC  is  operated 
with  a  constant  steady- state  utilization  factor  by  controlling  the 
natural  gas  flow  to  the  stack  as  [8] 


where  Ufs  is  the  desired  utilization  in  steady-state.  Furthermore, 
because  the  fuel  processor  is  specially  considered,  the  relation¬ 
ship  between  a  small  change  of  stack  current  A /  and  a  small 
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change  of  hydrogen  input  in  A q1^  fed  to  the  SOFC  stack  can  be 
derived  as  [8]: 


a7h2 


N  2Kr 

--  - A  7=  A I 

2Fufs(  \  +  tf.y)  Mfs(l  +  tfs) 


(27) 


FromEq.  (4),  the  small-signal  relationship  between  Awf,  A 7, 
and  Agj?  about  their  initial  steady-state  values  Ufs,  Io,  and  0 
can  be  written  as 


*,  ^h2,o  .  ,  Mfs  .  in 

A/  =  A«f  +  2^;  A<?h2  (28) 

Under  constant  Wf  control,  substituting  Eq.  (27)  into  Eq.  (28), 
and  Eq.  (28)  can  be  rewritten  as 


Awf  = - A I 

Io  1  +  Tf.y 


(29) 


Fig.  4.  Cell  terminal  voltage  under  various  current  and  hydrogen  flow. 


where  Iq  the  initial  stack  current,  is  related  to  the  initial  0 
through  Eq.  (4). 


5.  Simulation 


by 


RMSE  = 


(y(k)  -  y(k))2 


(30) 


In  this  section,  we  present  numerical  experiments  to  show 
the  validation  of  the  proposed  nonlinear  predictive  control 
scheme  based  on  the  GA-RBF  neural  network  model  of  the 
SOFC. 

For  the  purpose  of  identification,  the  input  signals  of  the 
dynamic  physical  model  developed  in  Section  3  are  uniformly 
random,  including  the  hydrogen  flow  rate  (0-1.2  mol  s-1)  and 
the  current  (0-800  A).  To  obtain  values  at  integer  time  points,  the 
fourth-order  Runge-Kutta  method  was  used  to  find  the  numer¬ 
ical  solution  to  the  dynamic  physical  model  in  the  simulation. 
A  set  of  3000  data  was  collected  from  the  simulation.  The  first 
1000  data  were  used  for  the  identification  of  GA-RBF  model, 
while  the  remaining  2000  data  were  used  for  validation  pur¬ 
poses. 

In  order  to  reduce  the  numbers  of  the  parameters  and  improve 
the  speed  of  program  debug,  the  hidden  layer  of  the  RBF  neural 
network  is  chosen  2  nodes.  After  many  trials,  a  population  size  of 
50,  a  crossover  probability  of  0.4,  and  a  mutation  probability  of 
0.001  are  used.  The  optimized  initial  values  of  the  parameters  are 
shown  in  Table  2.  After  the  optimized  initial  values  of  the  three 
parameters  are  obtained,  we  utilize  the  gradient  descent  learning 
algorithms  to  adjust  them.  Root  mean  square  error  (RMSE)  is 
employed  here  to  evaluate  modeling  results,  which  is  calculated 

Table  2 

The  optimized  initial  values  of  widths,  centers  and  output  weights 

bi  0.4856 

b2  1.1062 

cii  0.5678 

ci2  -0.3561 

c2i  —0.8409 

C22  1.032 

Wi  0.4567 

w2  -0.2365 


where  A  is  the  number  of  sample  data  from  the  dynamic  physical 
model  of  SOFC,  y(k )  is  the  predictive  output  of  GA-RBF  model, 
y(k)  is  the  output  of  the  dynamic  physical  model.  Under  various 
hydrogen  flows  and  current,  the  voltage  identification  model  is 
shown  in  Fig.  4.  The  RMSE  of  output  voltage  obtained  is  1 .2584 
and  1.1836,  respectively.  The  result  shows  that  the  GA-RBF 
neural  network  can  approximate  the  behavior  of  the  physical 
SOFC  model  with  good  accuracy..  The  control  input  is  shown 
in  Fig.  7 

In  normal  working  condition,  the  current  of  the  SOFC  system 
is  300  A.  The  steady  output  of  the  voltage  is  330  V.  Assuming 
at  t  -  200  s,  a  load  disturbance  causes  the  stack  current  to  have  a 
step  change  (from  300  to  255  A),  and  at  t= 550  s,  a  load  distur¬ 
bance  causes  the  stack  current  to  have  a  step  change  (from  255 
to  345  A).  The  series  of  step  in  stack  current  input  are  shown 
in  Fig.  5.  The  predictive  controller  is  used  to  adjust  the  voltage 
to  its  steady  value.  In  this  paper,  let  the  predictive  horizon  be 

350 1 - « - , - , - , - , - , - , - * - , - 

340  - 
330- 
320- 
g  310- 

|  300  - 

3  290  - 
280  - 
270  - 
260  - 

250 1 - 1 - l - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

0  100  200  300  400  500  600  700  800  900  1000 

Time  (s) 


Fig.  5.  Step  increase  in  cell  current. 
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Fig.  6.  Trajectories  of  SOFC  voltage  by  MPC. 


Fig.  9.  Fuel  utilization  response  by  constant  fuel  utilization  control. 


m=  10,  and  the  control  horizon  n  =  4.  We  can  get  the  tracking 
curve  of  the  controlled  voltage  of  the  SOFC  system,  which  is 
shown  in  Fig.  6.  From  Fig.  6,  we  can  get  when  the  stack  cur¬ 
rent  fluctuates  the  voltage  is  changed  rapidly  and  then  returns 
to  the  reference  value.  Based  on  the  results,  we  can  get  the  fuel 


Fig.  7.  Curves  of  manipulated  variable  e/)"  . 
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Fig.  8.  Fuel  utilization  response  by  MPC. 
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Fig.  10.  Curves  of  manipulated  variable  qf. 


utilization  response  curve  of  the  SOFC  system,  which  is  shown 
in  Fig.  8.  The  result  shows  that  the  predictive  control  scheme 
can  guarantee  the  fuel  utilization  of  the  SOFC  operating  within 
a  safe  range  when  load  changes. 

Now  we  use  the  same  load  disturbances  as  before  and  use 
the  constant  fuel  utilization  control  method.  The  tracking  curve 
of  the  controlled  fuel  utilization  of  the  SOFC  system  is  shown 
in  Fig.  9.  And  the  manipulated  variable  qf  is  shown  in  Fig.  10. 
Comparing  Fig.  8  with  Fig.  9,  we  can  get  the  conclusions:  the 
constant  fuel  utilization  control  method  is  simpler;  however  the 
fuel  utilization  uf  has  a  larger  excursion  from  the  safe  range 
when  the  stack  current  changes  from  255  to  345  A  at  7  =  550  s. 

6.  Conclusions 

Fuel  utilization  is  one  of  the  most  important  controlled 
variables  in  the  SOFC  system.  In  order  to  guarantee  the  fuel 
utilization  would  operate  within  a  safe  range,  a  nonlinear  pre¬ 
dictive  control  algorithm  based  on  a  GA-RBF  identification 
model  is  proposed  to  control  the  SOFC  output  voltage  as  a 
proper  constant.  By  comparing  the  constant  fuel  utilization 
control  method  with  the  nonlinear  model  predictive  controller 
scheme,  the  simulation  results  show  that  the  MPC  is  com- 
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plex,  however,  it  guarantees  the  fuel  utilization  performs  more 
safely. 
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